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Abstract 

We derive new bounds on decaying Dark Matter from the gamma ray measure- 
ments of (i) the isotropic residual (extragalactic) background by Fermi and (ii) the 
Fornax galaxy cluster by H.E.S.S.. We find that those from (i) are among the most 
stringent constraints currently available, for a large range of DM masses and a variety 
of decay modes, excluding half-lives up to ~ 10 26 to few 10 seconds. In particular, 
they rule out the interpretation in terms of decaying DM of the spectral features in 
PAMELA, Fermi and H.E.S.S., unless very conservative choices are adopted. We also 
discuss future prospects for CTA bounds from Fornax which, contrary to the present 
H.E.S.S. constraints of (ii), may allow for an interesting improvement and may become 
better than those from the current or future extragalactic Fermi data. 
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1 Introduction 



The possibility that Dark Matter (DM), which constitutes most of the matter in the Uni- 
verse, consists of a particle that actually decays on a very long time scale has attracted 
much attention lately. This is because the decay time scale Td ec can be taken to be 'short' 
enough that the decay products give signals in current high energy cosmic ray experiments. 
Namely, if Td cc — few ■ 10 26 sec, decaying DM can be invoked to explain the excesses in 
the fluxes of positrons and electrons measured by PAMELA, Fermi and H.E.S.S. [1]. On the 
other hand, this value of T^ ec is so much longer than the age of the Universe that the slow 
decay does not make a dent in the overall cosmological DM abundance and does not spoil 
the agreement with a number of astrophysical and cosmological observations [2] x . 

Irrespectively of this recent activity spurred by the charged CR anomalies, decaying DM 
is an interesting subject by itself. A long time studied case is the one of supersymmetric 
gravitino DM [4], which is unstable due to R-parity violation. More generically, in several 
particle physics models, high-scale suppressed operators may naturally mediate the decay 
of DM (for a miscellaneous list of references see [5, 6, 7]). 

From the phenomenological point of view, the main feature of decaying DM with respect 
to the 'more traditional' annihilating DM is that it is less constrained by neutral messenger 
probes (essentially gamma rays, but also neutrinos) originating from dense DM concentra- 
tions such as the galactic center, the galactic halo or nearby galaxies. The reason is simple 
and well-known: while the signal originating from annihilating DM is proportional to the 
square of the DM density, for decaying DM the dependence is on the first power; as a conse- 
quence, dense DM concentrations shine above the astrophysical backgrounds if annihilation 
is at play, but remain comparatively dim if DM is decaying. Decaying DM 'wins' instead, 
generally speaking, when large volumes are considered. This is why in the following we will 
focus on targets as large as galaxy clusters or, essentially, the whole Universe. 

On the observational side, the Fermi and H.E.S.S. telescopes are making unprecedented 
progress in the field of gamma-ray astronomy, producing measurements of many different 
targets including those of interest for decaying DM. It is therefore a good time to assess 
the current status of the latter. To this aim, we will compute the predicted signal from 
decaying DM, for a variety of decay channels (limiting ourselves to 2-body channels, in 
order to remain as model-independent as possible), and compare it to the gamma-ray 
measurements, deriving constraints on the decay half-life. In particular, in this paper, we 
make use of two distinct probes: 

(I) The isotropic residual gamma ray flux recently measured by Fermi [8], which now 
extends from about 200 MeV up to 580 GeV. The high-energy portion of this mea- 
surement is labelled as 'preliminary' by the Fermi collaboration, so it should be used 
with care. However, we note that the data are based on solid premises, since they are 
obtained with the same procedure already published in [9] with an enlargement of the 
dataset. We will anyway show later the effect of considering or not the preliminary 
portion of the data. 

(II) The recent observation in gamma rays of the Fornax galaxy cluster by H.E.S.S. [10]. 

1 See instead [3] for cases in which the decay of DM may actually help cosmology or astrophysics, but 
for much shorter decay time scales, excluded by the bounds we will present below. 
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We choose these probes for a number of reasons, besides the obvious fact that they are 
among the most recent ones. For what concerns the isotropic flux (I), it is known since a long 
time that it represents a powerful testbed for decaying DM (see for instance [11, 12, 13]). For 
what concerns (II), we are again motivated by the fact that large virialized objects such 
as galaxy clusters are a promising target for decaying DM (as briefly mentioned above) 
and, more in particular, by the fact that indeed stringent constraints have been recently 
derived using Fermi data [14, 15]: on the basis of this we explore the constraining power 
of a complemetary observatory such as H.E.S.S.. This is particularly interesting because 
upcoming Cerenkov telescopes offer very promising prospects of improvement in the near 
future (as opposed to a space based gamma ray observatory such as Fermi, which will at 
most increase its statistic by a factor of order 2). Indeed, we will also study the sensitivity 
of CTA, the upcoming large Cerenkov Telescope Array [16]. 

Incremental improvements to the data and the tools relevant for our analysis also comes 
from other aspects: 

(a) New data on charged Cosmic Rays (CR), which have been recently presented. Fermi [17] 
has measured the e + /(e + + e~) fraction, confirming the notorious rise exposed by 
PAMELA in 2008 [18] and extending the measured spectrum to larger energies (about 
200 GeV); Fermi has also presented updated measurements of the total flux of (e + + 
e~) [19] as well as the measurement of the separated fluxes of e + and of e~ [17]; MAGIC 
has published [20] the measurement of the total flux of (e + + e~) too; finally, PAMELA 
has presented results on the flux of pure e~ [21], of p [22] and of p [23]. 

(b) Electro Weak (EW) corrections to the DM generated fluxes. As discussed in a number of 
recent works [24], the emission of the EW gauge bosons (W ± , Z) from the SM particles 
emerging from the DM decay process can significantly modify the phenomenology, (i) 
The corrections are particularly relevant for large DM masses (above a TeV); (ii) they 
can alter significantly the fluxes of charged cosmic rays and gamma rays, both in their 
spectral shape and in their amplitude, affecting especially the low energies portion [25]. 
We include EW corrections in all our computations. 

(c) An improved propagation scheme for e ± in the Galaxy. In semi-analytic treatments of 
charged CR propagation in the galactic environment, usually employed in DM stud- 
ies, the energy losses of e + and e~ are generally considered as space-indepedent and 
their behaviour with energy E is approximated as E 2 . However, a space dependence 
obviously arises from the different conditions (of magnetic field and of ambient light) 
in different points of the galaxy and an energy dependence different from E 2 occurs in 
the full Klein-Nishina regime of Inverse Compton scattering. We make use of the fluxes 
of e ± as provided in [26], which include these proper space- and energy-dependences. 
This impacts both the analysis of the fits (mildly) and the fluxes of ICS gamma rays 
from the galactic halo (more significantly). 

Some recent studies have performed work related to ours [27, 28, 29, 30, 31, 32, 14, 33, 
15]. The most closely related, as far as the isotropic gamma ray flux is concerned, is [29], 
in which a subset of us derived constraints on decaying DM with the data and the tools 
available at that time. We improve, with respect to that work, by including the updated 
data in I and by the points (a), (b) and (c) discussed above. There are also some other 
minor analysis differences on which we will comment in the text. Among other works, 
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ref. [15] made a particularly detailed analysis: we consider the more recent data in I and 
II above and we include the Inverse Compton Scattering contribution to the gamma ray 
flux in the galactic component of I; the points (a), (b) and (c) are new in our paper; we 
also consider a few more decay channels. For the cases which overlap with this and other 
studies, we will present a comparison in Sec. 5.1. 

The rest of this paper is organized as follows. In Sec. 2 we update charged CR fits by 
including (a), (b) and (c). In Sec. 3 we discuss the calculation of the constraints from the 
isotropic residual 7-ray flux, while in Sec. 4 we discuss those from the Fornax cluster. In 
Sec. 5 we present the combined results. In Sec. 6 we present our conclusions. 

2 Update of the decaying DM fits to charged CR 
anomalies 

As mentioned in Sec. 1, the anomalous PAMELA, Fermi and H.E.S.S. data in e + and (e + +e~) 
have been interpreted in terms of DM decay. We recall here briefly the main features of 
the experimental data and of their DM interpretations, without entering in the details of 
any specific particle physics model. 

We use the following data sets: 

o PAMELA positron fraction [18], selecting only points with E > 20 GeV, in order both 
to avoid the low energy region affected by the uncertainty of solar modulation and to 
have a consistent overlap with the Fermi positron fraction data points (the low energy 
PAMELA data have very small error bars that would overconstrain the fit). 

o Fermi positron fraction [17]. 

o Fermi (e + + e _ ) total flux [19], provided in the low energy (LE) and high energy (HE) 
samples. 

o H.E.S.S. (e + + e~) total flux [34, 35], also provided in a lower energy portion and a 
higher energy one. 

o MAGIC (e + + e~) total flux [20], which however consists of only 6 data points, with 
error bars larger than those of Fermi and H.E.S.S. at the same energy and therefore 
has no effect on the global fit. 

o PAMELA p flux [23]. 

We perform the fit to these data using the DM generated fluxes as provided in [26], which 
include the features discussed in points (b) and (c) above. In looking for the best fitting 
regions, we scan over the propagation parameters of charged cosmic rays and over the 
uncertainties on the slope and normalization of the astrophysical electron, positron and an- 
tiproton background. Note that a simple power law is not expected to necessarily provide a 
good background model, especially given that the high energy part of the electron spectrum 
should be dominated by the contribution from local sources. On the other hand, due to 
the limited number of points available at high-energy and their relatively large statistical 
error, the background power-law flux is significantly constrained by the lower energy data 
and the assumption is sufficient in providing a good fit, when combined with the bumpy 
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spectrum from dark matter. We have assumed a NFW profile for the galactic DM halo (the 
specific parameters of which are discussed in Sec. 3), but other choices would have given 
almost indistinguishable results. We refer to [29] and [36] for more details. 

Fig. 1 presents the spectra results for one specific case: the best fit point for the decay 
mode DM—)- / u + ^ . The allowed regions on the plane M-QM~ T deo for the same channel, are 
instead shown in fig. 2. We show the 95.45 % and 99.999 % C.L. regions for the fit to 
positron and antiproton data only (green and yellow bands) and for the whole datasets, i.e. 
including (e + + e~) (red and orange blobs). In fig. 3, which we will discuss in more detail 
later, we show the global allowed regions for all the channels that we consider. 

The typical decay time scales that are required for the global fit are of the order of 10 26 
to 10 27 seconds. Moreover, as well known [36], only 'leptophilic' channels allow a global fit: 
for the quark and gauge boson channels, the few TeV decaying DM needed by (e + + e~) is 
in conflict with p data. 

The impact of the new data and of the improved analysis tools on the identification of 
the best fit DM properties is overall not big. For instance, in fig. 2 we show the dotted 
contours of the formerly preferred regions (taken from [29]), which are not very different. 
However, some features can readily be identified: 

- The fit region for the positron and antiproton signals (yellow and green, see fig. 2) 
starts now at larger masses, due to the fact that Fermi e + data reach higher energies 
than the former ones from PAMELA. 

- An allowed fit region now appears also for the channel DM-) e + e~ (see fig. 3); it was 
not present in previous analyses such as [29]. This is due to a number of concurring 
reasons, among which the fact that the inclusion of the EW corrections (b) and the 
refined propagation (c) smoothens out the e + + e~ spectrum and allow a decent fit 
to the ~ 1 TeV hump, the fact that we now remove from the fit the low energy e + 
PAMELA datapoints and, also, that we are now employing for Fermi data the full error 
bars (slightly larger than before). 

- EW corrections imply the presence of a non-zero p flux even for leptonic decay chan- 
nels (see e.g. the third panel of fig. 1). However, the yield is not enough to appreciably 
affect the fit regions. 

A comment on the selection of the datasets and their compatibility is also in order. We use 
the full set of data listed at the beginning of this section, but we also checked that dropping 
the data by PAMELA on the positron fraction, i.e. using only data from Fermi and H.E.S.S. 
(on the positron fraction and the pure e + and e~ fluxes), does not significantly change the 
fit regions. Similarly, the regions are not significantly modified if we adopt for the positron 
astrophysical background the recent determination in [37] (instead of the one conventionally 
used, from [38, 39]). This is expected, since the predictions in [37] and in [38, 39] are very 
similar in the range of energies in which we are interested. 

3 Isotropic gamma ray flux 

The measurements in (I) by the Fermi satellite correspond to the (maximal) residual, 
isotropic gamma-ray flux present in their data. Its origin can be in a variety of differ- 
ent phenomena, both in the form of unresolved sources and in the form of truly diffuse 
processes (see [8] and reference therein). 
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Figure 1: An example of the signals in different channels from a decaying DM candidate 
that fits the charged CR anomalies: we show the positron fraction (upper left), all leptons 
(upper right), antiproton (lower left) and isotropic gamma ray (lower right) spectra for the 
best fit candidate DM — )■ , namely Mdm = 3760 GeV, t& cc = 1.6 10 26 sec (marked by a 
white cross in fig. s 2 and 3). The 'PAMELA 2008' datapoints are reported for completeness, 
but we use only the more recent 'PAMELA 2010' in the fits (see text for details). In each 
panel, the DM contributions are dashed and the astrophysical background is shaded gray. 
In the lower right panel, the dotted (dashed) lines refer to the flux neglecting (including) 
extragalactic absorption, see Sec. 3. 
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Figure 2: Illustrative example (for the channel DM-^ / u + 'fjT ) of the impact on the fit and 
constrained regions following from different assumptions and choices, as discussed in the 
text. In fig. 3, on the other hand, we will show only the final regions and our fiducial 
constraints. 



DM decays can also contribute to this isotropic flux, with two terms: 1) an extragalactic 
cosmological flux, due to the decays at all past redshifts; 2) the residual emission from the 
DM halo of our Galaxy. The former is of course truly isotropic, at least as long as one 
neglects possible nearby DM overdensities. The latter is not, but its minimum constitutes 
an irreducible contribution to the isotropic flux. In formulae, the predicted differential DM 
flux that we compare with Fermi isotropic diffuse 7-ray data is therefore given by 
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For typical DM decay channels and for any DM halo profile, we find that the two contribu- 
tions are of comparable amplitude, as it can be seen e.g. in the lower right panel of fig. 1. 
We now move to discuss the calculations of the two contributions separately, reproducing 
and updating the discussion of [29]. 

The extragalactic flux is given, in terms of the Earth- measured photon energy E y , by 
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where Tdec = t^} c is the decay rate. Here the Hubble function H{z) = if \AUf(l + z) 3 + ^a, 
where H Q is the present Hubble expansion rate. Qdm, and fl\ are respectively the dark 
matter, matter and cosmological constant energy density in units of the critical density, 
p C) o- The gamma ray spectrum dN/dE^, at any redshift z, is the sum of two components: 
(i) the high energy contribution due to the prompt 7-ray emission from DM decays and 
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(ii) the lower energy contribution due to Inverse Compton Scatterings (ICS) on CMB pho- 
tons of the e + and e~ from those same decays. Using eq. (2), the extragalactic flux can 
therefore be computed in terms of known quantities for any specified DM mass Mdm and 
decay channel. In particular, since the flux essentially consists of the integral over p, at the 
first power, the result is not affected by the formation of DM halos, the history and the 
properties of which are highly uncertain. This is in contrast to the case of annihilating DM, 
where instead these issues affect the predictions by orders of magnitude (see e.g. [40, 41]). 
We take the resulting fluxes from [26], but we improve the treatment of the effect of the 
finite optical depth of the Universe, as we now move to discuss. 

The factor e^ T<yE ^ ,z ^ in eq. (2) accounts for the absorption of high energy gamma rays 
due to scattering with the extragalactic UV background light: essentially, a high energy 
photon from DM hits a lower energy UV photon and produces an energetic electron (or 
positron) which, in turn, makes Inverse Compton Scattering on the CMB and therefore 
reinjects a continuum of lower energy 7-rays. The process amount to a subtraction of 7-ray 
flux from the higher energies and a redistribution towards the lower part of the spectrum. 
As an element of novelty in this analysis, we take this process fully into account, using the 
energy and redshift dependent optical depth of the Universe (denoted t(E 1 , z)) of [26]. The 
effect is sizable at energies E 1 > 100 GeV: it can reduce the high energy flux by about one 
order of magnitude. Indeed, in fig. 1, lower right panel, we show with dotted lines the fluxes 
computed neglecting absorption and with dashed lines the effect of including it. Summing 
over the galactic component (which is of course not affected by absorption) reduces the 
impact of the effect, which can however remain sizable. This will affect the constraints 
derived below, especially for channels which feature a large prompt contribution. 

The flux from the galactic halo, coming from a generic direction dfl, is given by the well 
known expression 



i.e. as the integral of the decaying DM density piling up along the line of sight individuated 
by the direction dfl. Here Phaio is the DM distribution in the Milky Way, for which we will 
always take a standard Navarro-Frenk- White [42] profile 



with parameters r s = 24.42 kpc and p s = 0.184 GeV/cm 3 [26]. In principle one could 
consider other choices of profiles and parameters, but, for the case of decaying DM and 
considering that we will not be interested in the regions of the Galactic Center (GC) where 
profiles differ most (quite the opposite: we will focus on the anti-GC), these choices make 
a negligible difference. The coordinate r, centered on the GC, reads r(s, ip) = (r| + s 2 — 
2r scos-ip) 1 ^ 2 , where r = 8.33 kpc is the most likely distance of the Sun to the GC and 
ip is the angle between the direction of observation in the sky and the GC. 

The spectrum dN/dE^ consists again of two components: the prompt one and the ICS 
one. We compute both using the tools in [26]. In particular, for the ICS flux we use 
the generalized halo functions for the IC radiative process provided there, which take into 
account the full energy losses discussed in (c). This is another refinement with respect 
to [29], which had computed the e energy losses towards the anti-GC including CMB 
only: it was thus missing synchrotron losses (which can account for up to roughly a third 




(3) 




(4) 
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of the total losses, in the range 8.33 kpc < r < 15 kpc) and ICS losses on local ambient 
light other than the CMB (which are instead less important). The treatment that we now 
adopt represents the most refined one we can afford, at the current state-of-the-art of these 
semi-analytical computations. 

As indicated in eq. (1) and in the discussion above it, we need to determine the minimum 
of the flux in eq. (3). For the prompt contribution, the minimum is obviously located where 
the line-of-sight density of DM is the lowest, just because the 7-ray emission density traces 
the DM density by definition. This coincides with the direction of the anti Galactic Center 
(anti-GC). For what concerns the ICS contribution, on the other hand, the situation is 
more complicated. In this case the source of 7-rays is the population of DM-originated 
e , and therefore the minimum depends on what is assumed for the propagation and final 
distribution of the latter ones in the galactic halo. For instance, assuming a thin diffusive 
halo with L = 1 kpc (see e.g. [26] for a discussion and an overview of standard values) implies 
that the layer of ICS emitting e ± is particularly shallow in the zenith (and nadir) direction 
and therefore the ICS 7-ray flux is minimal along lines of sight pointing directly above (and 
below) the solar system. On the other hand, for a thick L = 15 kpc diffusive halo, the 
ICS emission from the directions orthogonal to the galactic plane is greatly enhanced and 
the minimum is located again at the anti-GC. In view of this complicated morphology, we 
choose to locate the minimum always at the anti-GC: 



This is supported, along the lines of the discussion above, by a number of considera- 
tions, among which: (i) thin diffusive halos are anyhow disfavored by other arguments 
(see e.g. [43]); (ii) a number of other uncertain astrophysical variables would enter in a de- 
tailed determination of a true minimum: for instance if the radial size of the diffusive halo is 
smaller than usually assumed (~20 kpc), then the ICS flux at the anti-GC is reduced; (iii) 
this discussion only affects the galactic ICS contribution (not the galactic prompt and the 
extragalactic contributions), thus limiting the impact on the derived total flux and having 
an importance (for the bounds we will derive below) only for those configurations in which 
the ICS saturates the constraints. 

3.1 Deriving the constraints 

We thus compute G?$ isotrop i C /(i.E' 7 as discussed above and we compare it with the Fermi data 
of (I). We show an example of such a comparison in Fig. 1, lower right panel. It is pretty 
clear (even by inspection, in the case of this example) that the DM signal does not agree 
in shape with the data, which are instead well fit by a simple power law [8]. For this 
reason, we are driven to derive constraints on the maximum DM signal, and therefore the 
minimum Td ec , admitted by the data. There are however several possible ways to compute 
such constraints. We discuss a few of them in the following and we illustrate them in the 
specific example of fig. 2. 

o DM signal only. The most conservative option consists in demanding that the DM 
signal alone (eq. (1)) does not exceed any single one of the Fermi data points by 
more than a chosen significance, which we set at 3a. This is overly conservative for 
a number of reasons. First of all, since it supposes a vanishingly small astrophysical 
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background. A simple inspection of the DM prediction proves that this assumption is 
physically untenable; the bumpy shape of the signal is so different from the featurless 
observations that the (bulk of the) latter must be explained by astrophysical processes, 
leaving only a subleading role for signals from DM. Second, given the smooth nature 
of the DM signal and of the data, it is clear that an excess of 3a in one point is often 
accompanyed by similar excesses in neighboring points, so that the 'global significance' 
of the exclusion is actually higher. For instance, the best-fit model example in fig. 1 
is only barely excluded by such a procedure, despite the clear tension with several 
data points. 

o DM signal + power-law background. A more realistic option consists in assuming 
that the astrophysical background consists of a power-law and demand that the sum 
of astrophysical background and DM signal does not exceed a chosen global signifi- 
cance. Note that one expects several potential sources to contribute to this flux, like 
unresolved blazars, star-forming galaxies or electromagnetic cascades from ultra-high 
energy cosmic ray losses. In general, a combination of them is not expected to pro- 
duce (a priori) an exact power-law, which may be inaccurate also for some of the 
different contributions taken alone. Nonetheless, within the current precision we find 
that the data of E^d§/dE 1 dVl alone are very well described by a power law with 
index -0.41 (very similar to the one found in resolved blazars) and normalization 1.02 
10~ 5 GeV/(cm 2 sec sr). Next, for a given Mom we add the DM signal, whose normal- 
ization is controlled by rd ec . We let the normalization of the power law background 
vary within a factor of 2 from the central value specified above and its index within 
and -1 (this choice is astrophysically plausible, although varying the parameters in a 
broader range would not change the results shown). We then compute the x 2 to the 
data, marginalizing over these parameters. We compute in this way 95% CL. limits 

on r dec . 

o DM signal + unconstrained background. A variation of the procedure described in the 
previous point consists in avoiding to commit on the functional form of the background 
and assume that it consists of an arbitrary function which describes well the data 
(namely: it is consistent with each data point within a small significance, e.g. la). 
We then demand that, adding the signal, the A% 2 does not exceed a given value, 
which in this case we fix very conservatively at 25. For the practical case at hand, 
however, this procedure and the previous one yield similar results, as one can see on 
the example in fig. 2. This is not surprising, since the isotropic residual datapoints 
are indeed very well described by a power law. 

The procedure 'DM signal + power-law background' is the one that we will adopt to obtain 
the fiducial constraints shown in our results (see fig. 3). This procedure is 'fiducial' also 
in the sense that it matches the analysis we do for charged CR anomalies (see Sec. 2) and 
therefore fit regions and constraints are essentially consistent with each other in fig. 3. 

Before moving on, we discuss how the constraints would be modified if we removed from 
the analysis the datapoints above 100 GeV, that are labelled as 'preliminary' by the Fermi 
collaboration, i.e. if we limited the data to those published during 2010 in [9] only. In fig. 2 
the dotted gray line corresponds to how the fiducial constraints would be modified. We see 
that, as expected, the limits become looser by a factor of a few. 
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4 Fornax cluster gamma ray flux 



Galaxy clusters are the largest gravitationally bound structures in the universe, 80% of 
their total mass being in the form of Dark Matter. Although they are located at much 
larger distances than other popular targets such as dwarf satellites of the Milky Way, they 
turn out to be attractive environments to search for DM due to their predicted high DM 
luminosity. In practice, which cluster is the most promising depends on a number of factors, 
and most notably the astrophysical background. Indeed, standard astrophysical gamma- 
ray emission is expected, both from high energy cosmic rays interacting in the intra cluster 
medium (the dominant contribution) or from electrons in the ambient radiation field (a 
subdominant contribution, see [44] for a review). Despite these predictions, no gamma-ray 
from these processes has been detected so far [45, 46, 47, 48, 49, 10]. Radio galaxies lying 
at the gravity center of clusters such as M87 in Virgo and NGC1275 in Perseus are also 
7-ray emitters: these have actually been directly observed. 

A few galaxy clusters have been observed by H.E.S.S.: the most attractive of them for 
DM searches are Virgo, Coma and Fornax. Although Virgo is closeby, the high energy 
gamma-ray emission from M87 prevents searches in the inner region. Coma is among the 
most massive galaxy clusters but it has been shown that it is not a privileged environment 
because of a relatively high CR-induced gamma-ray emission [52, 53]. Fornax is located at 
19 Mpc [50] and shows favorably-low expected astrophysical background (see, for instance, 
[51, 52, 53]). Given its location near the tropic of Capricorn, the H.E.S.S. instrument is best- 
suited to observe Fornax with respect to other currently-operating Imaging Atmospheric 
Cerenkov Telescopes (IACTs). 

The predicted DM 7-ray flux from Fornax can be easily obtained by integrating Eq.(3) 
over the observational solid angle Afl, and of course replacing p halo with ^Fornax, the DM 
distribution in the cluster. 

There are a few different methods to determine pFomax (and in general the DM content 
of galaxy clusters). A widely used approach is based on X-ray measurements on the grav- 
itationally bound hot intracluster gas. Assuming a NFW profile [42], the latter is entirely 
defined from the virial mass extracted form the HIFLUGCS catalog catalog [50] and the 
virial mass and the concentration relationship found in Ref. [54], in ACDM cosmology. 
Other choices for the profile are possible, e.g. the cored Burkert profile [55], which may 
be particularly motivated by the fact that baryon physics in the inner part of clusters may 
significantly alter the predictions of ACDM simulations: several processes invoking baryons 
in galaxy formation such as dynamical friction, AGN feedback or gas outflows, may flatten 
the DM cusps in DM cores [56, 57]. For the case of Fornax, various dark matter halo models 
have been considered in [10] (and references therein). However, as opposed to the annihila- 
tion case, the gamma-ray flux from decaying DM is less dependent on the DM distribution 
inside the target due to the simple linear dependence of the astrophysical factor with DM 
density. Varying over the profiles discussed in [10], the difference in the flux factor is less 
than a factor of 3 for a given opening integration angle. Another approach to determine the 
DM content of clusters consists in using dynamical tracers (see again [10] and references 
therein). 

In the following, we take as 'fiducial' the X-ray based determination of the DM profile, 
hereafter referred to as the RB02 profile [10], which also well agrees with tracer dynamics 
at large distances. In this case, the DM distribution is inferred from a X-ray concentration- 
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virial mass relationship 2 and the uncertainty on the mass determination is about 10% [54]. 
As an alternative choice, we also consider the DW01 profile of [10], which is fully based on 
the dynamical tracers method. This generates a predicted DM 7-ray flux ~3 times smaller 
than RB02. 

Next, we discuss the choice of AQ. In contrast to annihilating dark matter for which 
most often the smallest opening angle provide the most sensitive searches, at least for 
cuspy profiles, decaying dark matter searches require optimization of the opening angle 
to guarantee the highest signal-to- noise ratio. As it is straightforward from Eq.(3), the 
luminosity scales with the size of the solid integration angle. On the other hand, background 
is increasing as well. We find that the optimization of the signal-to-noise ratio versus the 
opening integration angle for the dark matter halo profile RB02 implies the best integration 
region to be 0.5°, which corresponds to a solid angle of AQ = 2.4 x 10~ 4 sr. 

We also mention that an additional component to the decaying Dark Matter spectrum 
needs in principle to be considered, particularly for leptonic final states: that of Inverse 
Compton emission (indeed, as opposed to dwarf galaxies, in clusters the electrons lose energy 
primarily through ICS on the ambient radiation field and produce additional gamma rays 
in the final state). We include it in the computations for the DM — > fi + fi~ channel, but, 
due to the energy working range of H.E.S.S., this component becomes important only for 
very large DM masses (above ~ 30 TeV) and its effect is therefore not visible in our results. 

4.1 Deriving the constraints 

In order to extract exclusion limits on the Dark Matter lifetime from gamma-ray astro- 
nomical observation with IACTs, the background needs to be determined to constrain the 
decaying DM luminosity. The background is calculated in a region referred to as the OFF 
region, and the signal region as the ON region. Both ON and OFF regions depend on the 
observation mode and are specific to the IACT instrument. As for the background level, 
it is taken simultaneously, i.e. in the same data-taking observing conditions, to the signal 
events in order to allow for the most accurate estimate. 

In the 'standard' ON-OFF method, one has to take into account the background flux 
in the astrophysical factor and subsequently in the upper limit calculation. For the ON 
region of 0.5° and considering a OFF region in an annulus around the ON region with an 
inner and outer radii of 0.5° and 1.0° respectively, this reduces the decaying DM flux in 
the ON region of about 40%. An 'improved' background estimate procedure called the 
template method allows to avoid this reduction. In such method the background events are 
determined in the ON region but from selecting hadron-like events (see [10] for details on 
the analysis procedure). We indeed make use of it in this work. 

H.E.S.S. has observed Fornax for a total of 14.5 hours at low zenith angle to allow for 
best sensitivity to low dark matter masses. 

2 The virial galaxy cluster mass Ma is usually defined as the enclosed mass within a radius where the 
density reaches A x p c , with p c the critical density of the universe and A taking values between 100 and 
500. Although M500 can be used to derive specific galaxy cluster properties, the X-ray concentration-mass 
relation is particularly well verified for A = 200 [58] . 



11 



5 Results and discussion 



Figure 3 presents the exclusion plots and therefore summarizes our main results. One can 
see that the constraints from the Fermi isotropic 7-ray data exclude decaying DM with a 
lifetime shorter than 10 26 to few xlO 27 seconds, depending on its mass and the precise 
channel. Therefore, in particular, they rule out the charged CR fit regions, for all the 
channels. As illustrated in the example in fig. 2, and inspecting fig. 1, adopting the more 
conservative constraint procedure may marginally reallow a portion of the fit regions, for the 
DM — » yU + /i~, but leaving a clear tension. On the other hand, removing the 'preliminary' 
data of [8] and keeping only those published in 2010 in [9] still allows to exclude the CR 
fit regions (as illustrated in fig. 2 for the the other channels are less critical and 

remain safely excluded). The constraints from Fermi rise gently as a function of the mass, 
essentially as a consequence of the fact that the measured flux rapidly decreases with energy. 
They also depend (mildly, a factor of a few at most) on the decay channel, as a consequence 
of the different 7-ray yield and the different shape of the DM signal in each channel. 

The constraints from H.E.S.S. Fornax remain subdominant, roughly one order of mag- 
nitude below the Fermi ones. However, for the case of the DM — > t + t~ channel, the bound 
also reaches the CR fit region and essentially confirms the exclusion. The constraints are 
cut at low masses by the dying out of the sensitivity in H.E.S.S. to low energy photons. 
They do not look competitive with respect to Fermi even for larger masses. 

In this study we have focussed on 2-body particle-antiparticle decay modes (DM — > 
£ + £~ ,qq,W + W~ , where I is a charged lepton), typical of scalar DM. We do not address 
3-body decays such as DM — > £ + £~ / y, since those are model dependent. For fermionic DM, 
decay channels such as DM — > £^W T are possible: these (in first approximation) can be 
analysed in our framework with a trivial combination of the DM — > £ + £~ and DM — > W + W~ 
channels. Another possibility which has been recently considered is mixed ('democratic') 
leptonic channels, such as DM — > 33% e + e~ + 33% + 33% t + t~: these can of course 

be approximately derived on the basis of our individual bounds. 

5.1 Comparison with existing bounds 

We here comment on the relative strength of the constraints in fig. 3 with respect to bounds 
from other analyses or other targets. 

With respect to the isotropic 7-ray constraints of [29], the bounds derived here are 
stronger by a factor of 2 to 3. The reasons of this are discussed at length in the text and 
essentially amount to: updated datasets, more refined DM analyses and the adoption of a 
more realistic constraint procedure (some of these effects pull towards a weakening while 
others for a strengthening). 

With respect to the work in [15], our constraints from the Fermi isotropic background 
are somewhat stronger than their corresponding ones (up to about a factor of 5 for the 
DM — > fi + fi~ channel), likely due to the several differences in the analysis discussed in 
Sec. 1. In addition, [15] presents bounds from the observation of the the Fornax cluster by 
Fermi: these are less stringent than our Fermi isotropic background constraints but more 
powerful than our H.E.S.S. -based constraints at moderate masses. At the largest masses, 
our H.E.S.S. -based constraints pick up and match theirs (for the bb channel), as expected 
from the different operational energy range of Fermi versus H.E.S.S.. The bounds from 
clusters other than Fornax are less powerful, according to the analysis in [15]. Ref. [59] 
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Figure 3: The regions on the parameter space MDM~ r dec that are excluded by the Fermi 
and H.E.S.S. constraints and that can be explored by CTA ; together with the regions of the 
global fit to the charged CR data, for dijjererfy decay channels. 



has derived constraints from several clusters using Fermi data too, but their procedure is 
questioned in [15], and in any case it is now superseded. On the other hand, the preliminary 
constraints shown (for the bb channel only) in [14], obtained with a combination of several 
clusters in Fermi, exceed our bounds by a factor of 2. 

Bounds from probes other than the isotropic flux and clusters do not generally achieve 
the same constraining power. E.g. recently [60] finds less stringent limits for the leptophilic 
channels, focussing on the Milky Way halo. 

Within the context of observations performed by IACT, we note that the decay lifetime 
constraints obtained with galaxy clusters are stronger than those from dwarf galaxies. Even 
for the ultra-faint dwark galaxy Segue 1, which is believed to be the most promising dwarf 
in the northern hemisphere 3 , the constraints are 2 orders of magnitude higher for a dark 
matter particle mass of 1 TeV [61]. 

We also estimate that the constraints that we derive are stronger than those that can 
come from neutrino observation of the Galactic Center (see for instance [62]). Precisely 
computing those ones for all channels and in the same analysis framework of this study, 
however, is beyond the scope of this work, and would also probably benefit from a better 
knowledge of the Icecube data (especially the reconstructed neutrino energy), which is not 
currently available. 

In summary: with the possible exception of the preliminary bounds from a combination 
of clusters by Fermi for the bb channel, the constraints that we derive from the isotropic 
7-ray flux are the most stringent to date. 

5.2 Prospects for improvements 

As we have seen, currently the constraints on decaying DM from the Fermi satellite are 
dominant with respect to those from H.E.S.S.. However, while the former may increase its 
statistics by at most a factor of a few, for the latter there are prospects of developments in 
the mid-term future. 

The next-generation IACT will be a large array composed of a few tens to a hundred 
telescopes located on two sites, one in each hemisphere [63]. The goal is to improve the 
overall performances of the present generation: one order-of-magnitude increase in sensi- 
tivity and enlarge the accessible energy range both towards the lower and higher energies 
allowing for an energy threshold down to a few tens of GeV. From the actual design study 
of CTA, the effective area will be increased by a factor 10 and a factor 2 better in the hadron 
rejection is expected. 

The calculation of CTA sensitivity is performed with a background-only hypothesis, 
where the non-detection of a signal by CTA would imply that the signal in the ON region 
consists of only misidentified hadron showers. The background level is estimated by in- 
tegrating Eq. (4) of [64] after multiplying it by the CTA effective area. The upper limit 
on the number of gamma-ray events is calculated at 95% C.L. according to the method 
from [65] assuming five OFF regions. This calculation follows the methodology described 
in [64], from which the CTA effective area is extracted. The 95% C.L. sensitivity on the 

3 Ultra-faint dwarf galaxies are particularly interesting due to their high mass-to-light ratio. However 
the nature of some of them is still under debate due to their similarity of their properties with globular 
clusters. Due to their low surface brightness and a few tens of member stars, their dark matter content is 
still subject to large uncertainties. 
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decay lifetime is then given by 

r 95%C.L. = 47T M DM A^.L. 

La dn Ls ds P**Ur{*, n Tobs J Mdm/2 A cta (£ 7 ) §-{131 dE, ' 

where ]\[^ 5%C - L - is the limit on the number of gamma-ray events, Acta is the CTA effective 
area and T obs the observation time. 

Figure 3 shows the 95% C.L. sensitivity of CTA on the decay lifetime for the RB02 halo 
profile for 50h observation time and AQ = 2.4 x 10~ 4 sr. 

To conclude, we also mention that a technique which could allow for significant im- 
provements in the exploration of the parameter space of decaying DM using clusters is the 
one of stacking the observation of a large number of different clusters, as recently discussed 
in [66]. The authors find that improvements of up to 100 can be theoretically achieved, 
albeit this factor is ~5 for more realistic background-limited instruments. 



6 Conclusions 

Decaying Dark Matter has come to the front stage recently as an explanation, alternative to 
annihilating DM, for the anomalies in CR cosmic rays in PAMELA, Fermi and H.E.S.S.. But, 
more generally, decaying DM is a viable possibility that is or can naturally be embedded 
in many DM models. It is therefore interesting to explore its parameter space in the light 
of the recent observational results. 

We discussed the constraints which originate from the measurement of the isotropic 
7-ray background by Fermi and of the Fornax cluster by H.E.S.S., for a number of decaying 
channels and over a range of DM masses from 100 GeV to 30 TeV. We improved the analysis 
over previous work by using more recent data and updated computational tools. 

We found (see fig. 3) that the constraints by Fermi rule out decaying half-lives of the 
order of 10 26 to few 10 27 seconds. These therefore exclude the decaying DM interpretation of 
the charged CR anomalies, (at least) for all 2-body channels, at least adopting our fiducial 
constraint procedure. The constraints by H.E.S.S. are generally subdominant. For the DM 
— > r + r~ channel, they can however also probe the CR fit regions and essentially confirm 
the exclusion. 

With one possible exception for the DM — > bb channel, the constraints that we derive 
from the isotropic 7-ray flux are the most stringent to date. 

We also discussed the prospects for the future Cerenkov telescope CTA, which will be 
able to probe an even larger portion of the parameter space. 
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